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1. Introduction 

Stochastic control is the study of dynamical systems subject to random per- 
turbations and which can be controlled in order to optimize some performance 
criterion. Historically handled with Bellman's and Pontryagin's optimality prin- 
ciples, the research on control theory considerably developed over these last 
years, inspired in particular by problems emerging from mathematical finance. 
The global approach for studying stochastic control problems by the Bellman 
dynamic programming principle has now its suitable framework with viscosity 
solutions concept : this allows to go beyond the classical verification Bellman 
approach for studying degenerate singular control problems arising typically in 
finance. On the other hand, the stochastic maximum principle finds a modern 
presentation with backward stochastic differential equations (BSDE) , which it- 
self led to a very important strand of research. Viscosity solutions and BSDE 
have a strong relation through their PDE representation, and stochastic control 
theory is a place where probabilists and PDE mathematicians meet together. In 
this survey paper, we give an overview of some of these modern developments 
of stochastic control by focusing mainly on controlled diffusion processes. Our 
chief concern is to derive some tractable characterization of the value function 
and optimal control. We do not discuss the theoretical existence problem of 
optimal controls, which was largely studied in the literature, and we refer to 
for a very recent overview on this subject. We include examples of typical 
applications with detailed explicit solutions. 

The paper is organized as follows. The next section formulates the problem, 
discuss the classical approaches and their limitations. We also present some 
other stochastic control problems of interest and their possible developments. 
In Section |3| we present the Bellman dynamic programming approach with vis- 
cosity solutions theory. Some original proofs are exposed and we indicate for 
some applications in finance how examples may be explicitly solved. Section 0] 
describes the modern exposition of stochastic control and maximum principle 
by means of BSDE. An application to the stochastic linear quadratic control 
problem and mean- variance hedging is explicited. In both sections, we give the 
key results on the characterization of optimality. In Section we discuss nu- 
merical issues and conclude in Section with some possible developments and 
open questions. 

2. The problem and discussion on methodology 
2.1. Problem formulation 

A standard control diffusion problem on finite horizon is formulated as follows. 
Let (f2, J 7 , P) be a probability space, T > a finite time, F = (J 7 t)o<t<T a filtra- 
tion satisfying the usual conditions, and W a d-dimensional Brownian motion 
defined on the filtered probability space (fi, F, P). We consider the controlled 
state process X valued in R™ and satisfying : 



dX s = b(s,X S7 a s )ds + a(s,X S7 a s )dW s . 



(2.1) 
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The control a = (a s )o<s<T is a progressively measurable process valued in the 
control set A, a subset of 1™. The Borelian functions b, a on [0,T] xl"xi 
satisfy the usual conditions in order to ensure the existence of a strong solution 
to 12.1fl . This is typically satisfied when b and a satisfy a Lipschitz condition on 
(i, x) uniformly in a, and a satisfies a square integrability condition. We denote 
by A the set of control processes a. Given (t,x) G [0,T] x R n and a € A, we 
then denote by {Xl' x ,t < s < T}, the unique strong solution to (|2.1|) starting 
from x at time s = t. As usual, we omit the dependance of X on a to lighten 
notations. 

We are then given two Borelian real-valued functions / and g respectively 
defined on [0, T] x M" x A and R™ and we define the gain function : 



J{t, x, a) 



f{s 1 X t ^,a s )ds + g{X t ^) 



for all (t, x) G [0, T] x K™ and a G A. Of course, we have to impose integrability 
conditions on / and g in order for the above expectation to be well-defined, e.g. 
a lower boundedness or linear growth condition. Now, since our objective is to 
maximize this gain function, we introduce the so-called value function : 



v(t, x) 



sup J(i, x, a) 



(2.2) 



For an initial state (t, x), we say that a G A is an optimal control if v(t,x) — 
J(t, x, a). 

Remark 2.1 We focus mainly in this survey paper on finite horizon control 
problems. The infinite horizon version of control problem (|2.2(l is formulated as 



sup E 

aeA 



e~ pt f{X s ,a s )ds 



where p > is a positive discount factor, enough, to ensure finiteness of the 
associated value function. Such problems are studied similarly by Bellman's 
optimality principle as for the finite horizon case. There is also no additional 
difficulties for considering more general discount factors : p — p(t, x, a), both in 
the finite or infinite horizon case. The formulation (|2.2(l captures some funda- 
mental structure and properties of optimal stochastic control, but there are of 
course many other types of control problems that are also important both from 
a theoretical and applied viewpoint. We shall list some of them and present 
possible developments later. 



2.2. Bellman's optimality principle 

Bellman's optimality principle, initiated by Bellman |S] and also called the dy- 
namic programming principle (DPP), is a fundamental principle in control the- 
ory : it formally means that if one has followed an optimal control decision until 
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some arbitrary observation time, say 9, then, given this information, it remains 
optimal to use it after 8. In the context of controlled diffusion described above, 
the DPP is mathematically stated as follows : 

(DPI) For all a € A and 9 £ % } t, set of stopping times valued in [t, T] : 



v(t,x) > E 



f{s,Xl' x , as )ds + v{6,Xl> x ) 



(DP2) For all e > 0, there exists a e £ A s.t. for all 9 e 7^ T : 



v(t, x) — e < E 



f( S ,X^,al)d S + v(e,X^ x ) 



(2.3) 



(2.4) 



Notice that this is a stronger version than the traditional version of dynamic 
programming principle, which is written as : for all 9 € T t .T, 



(DPP) 



v(t, x) 



sup E 



' ,0 

j f(s,Xl<*,a s )ds + v(e,Xt' x ) 



Although the DPP has a clear intuitive meaning, its rigorous proof is technical 
and has been studied by several authors and by different methods. We mention 
among them [HHj, [HI]; US HU or [HH]- However, it is rather difficult to 
find a precise reference with a self-contained and complete proof of the DPP 
in the above form DPI and DP2. Formal arguments are usually given but the 
technical part is often omitted or does not cover exactly the considered model. 
Actually, it is frequent to consider originally controlled diffusion model with 
strong solutions of SDEs, as formulated above. Then, in order to apply the dy- 
namic programming technique, one needs to consider a family of optimal control 
problems with different initial times and states along a given state trajectory. 
Hence, if AT is a state trajectory starting from x at time in a given filtered 
probability space (CI, T, F, P), then X t is a random variable, for any t > 0, in 
the original probability space (fi, T, P). However, since a control is F-adapted, 
the controller knows about all the relevant past information of the system up 
to time i, and in particular about X t . This means that X t is not uncertain for 
the controller at time t, or mathematically, almost surely deterministic under 
a different probability measure P(.\Tt). Thus, for using dynamic programming, 
we are naturally led to vary the probability spaces and so to consider the weak 
formulation of the stochastic control problem, for which one shows the dynamic 
programming principle. However, it is not clearly stated in the literature that 
the value functions of stochastic control problems under weak and strong formu- 
lations are equal. For this paper, we shall assume the validity of the DPP in the 
above form and mainly focus here on the implications for the characterization 
of the value function through the so-called Hamilton- Jacobi-Bellman (HJB) or 
dynamic programming equation. We end this discussion by mentioning that Bell- 
man's optimality principle goes beyond the framework of controlled diffusions or 
even Markov processes, and may also be applied for controlled semimartingales, 
see |5*T] . 
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2.2.1. The Hamilton- J acobi- Bellman (in)equation 

The HJB equation is the infinitesimal version of the dynamic programming 
principle : it describes the local behavior of the value function v(t, x) when 9 
is sent to t in DPP. It is formally derived by assuming that the value function 
is C 2 , applying Ito's formula to v(s,Xl' x ) between s = t and s = t + h, and 
then sending h to zero into DPP. The classical HJB equation associated to the 
stochastic control problem 11'. 21 is : 

dv 

-— (t,x) - sup [C a v(t,x) + f(t, x, a)} = 0, on [0,T)xR", (2.5) 

Ot a< =A 

where C a is the second-order infinitesimal generator associated to the diffusion 
X with constant control a : 

C a v — b(x, a). D x v + — tr (cr(x, a)<j' (x, a)D^v) . 

This partial differential equation (PDE) is often written also as : 

-—(t,x)-H(t,x,D x v(t,x),D 2 x v(t,x)) = 0, V(M)6[0,T)xM", (2.6) 

where for (t, x,p, M) € [0, T] x 1° x W 1 x S n (S n is the set of symmetric n x n 
matrices) : 



H(t,x,p,M) — sup 

aeA 



b(x, a).p + -tr (aa'{x, a)M) + f(t, x, a) 



(2.7) 



The function H is sometimes called Hamiltonian of the associated control prob- 
lem, and the PDE (|2.5|l or (|2.ti|l is the dynamic programming or HJB equation. 
There is also an a priori terminal condition : 

v(T,x) = g(x), VxeR™, (2.8) 

which results from the very definition of the value function v. We followed 
the usual PDE convention by writing a minus sign in front of the second order 
differential operator in the left side of l|2.5|) or (|2.6|) , which ensures that it satisfies 
the ellipticity condition, see (|3.3|l . 

Remark 2.2 The statement l|2.6|l of the HJB equation (in addition to the reg- 
ularity problem that will be handled with viscosity solutions) requires naturally 
the finiteness of the Hamiltonian H, which is typically satisfied when the set 
of controls A is bounded. Actually, when A is unbounded, it may happen that 
H(t, x,p, M) defined in (|2.7|l takes the value oo in some domain of (t, x,p, M). 
More precisely, assume there exists a continuous function G(t, x, p, M) on [0, T] x 
l"xl"x S n s.t. 

H(t,x,p,M) < oo G(t,x,p,M) > 0. 
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Then, from 1)2. we must have : 

G{t,x,D x v(t,x),D 2 x v(t,x)) > 0, (2.9) 
and - — (t,x) - H(t,x,D x v(t,x),D 2 x v(t,x)) > (2.10) 

Moreover, if inequality (|2.9|l is strict, then inequality (|2.10|l is an equality. Hence, 
we formally get the HJB variational inequality associated to the dynamic pro- 
gramming principle : 

Ov 

~~q1^^ x ) ~ H{t,x,D x v(t,x),Dlv(t,x)) , 

G(t,x,D x v(t,x),Dlv(t,x))] = 0. 

In this case, we say that the control problem is singular, in contrast with the 
regular case of HJB equation l|2.6|l . A typical example of a singular problem 
occurs when the control influences linearly the dynamics of the system and the 
gain function. To fix the idea, consider the one-dimensional case n = 1, A = 
R+, and 

b(x,a) = b(x) + a, <x(x,d) = a(x), f(t,x,a) — f(t,x) — Xa, 

for some A 6 R. Then, 

= {Kx) P+ \a{xfM + f{t,x) if"P+A>0 
v ' '- p ' ; 1 oo if -p + X < 0. 



The HJB variational inequality is then written as : 

flu 

= 0. 



--(t,x)-b(x)-(t,x)--a(x) —( t ,x),--(t,x) + X 

We shall give in the next section another example of singular control arising in 
finance and where the control is in the diffusion term. 

In a singular control problem, the value function is in general discontinuous 
in T so that (|2.8() is not the relevant terminal condition. We shall show how to 
derive rigorously the HJB equation (or variational inequality) with the concept 
of viscosity solution to handle the lack of a priori regularity of the value function, 
and also to determine the correct terminal condition. 



2.2.2. The classical verification approach 

The classical verification approach consists in finding a smooth solution to the 
HJB equation, and to check that this candidate, under suitable sufficient condi- 
tions, coincides with the value function. This result is usually called a verification 
theorem and provides as a byproduct an optimal control. It relies mainly on Ito's 
formula. The assertions of a verification theorem may slightly vary from prob- 
lem to problem, depending on the required sufficient technical conditions. These 
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conditions should actually be adapted to the context of the considered problem. 
In the above context, a verification theorem is roughly stated as follows : 

Verification theorem. Let w be a C 1 ' 2 function on [0, T) xl" and continuous 
in T, with suitable growth condition. Suppose that for all (t,x) & [0, T) x R™, 
there exists a(t, x) mesurable, valued in A s.t. w solves the HJB equation : 

sup [C a w(t, x) + f(t, x, a)] 

C &(t ' x) w{t,x) - f(t,x,a(t,x)) on [0,T)xM", 
together with the terminal condition 

w(T,.) = g on R n , 

and the S.D.E. : 

dX s = b{X s ,a{s,X s ))ds + a(X s ,a{s,X s ))dW s 

admits a unique solution, denoted X l / X , given an initial condition X t — x. Then, 
w = v and {a(s, Xl ,x ) t < s < T} is an optimal control for v(t, x). 

A proof of this verification theorem may be found in any textbook on stochas- 
tic control, see e.g. &\, |HI|> 0U, [SHI or EH- The first and most famous ap- 
plication in finance of this verification theorem for stochastic control problem 
is Merton's portfolio selection problem. This is the situation where an investor 
may decide at any time over a finite horizon T to invest a proportion a val- 
ued in A = R of his wealth X in a risky stock of constants rate of return /z 
and volatility a and the rest of proportion 1 — a in a bank account of constant 
interest r. His wealth controlled process is then governed by : 

dX s = X s (r + (n - r)a s ) ds + X s aa s dW s , 

and the objective of the investor is given by the value function : 

v(t,x) = m^E\U{X.y")] , (t,x) G [0,T] x R + , 

where U is a utility function, i.e. a concave and increasing function on R + . For 
the popular specific choice of the power utility function U(x) = x p ', withp < 1, it 
is possible to find an explicit (smooth) solution to the associated HJB equation 
with the terminal condition v(T, .) = U, namely : 

v(t,x) = exp (p(T - t)) x p , 

with p = 2a 2 T^p r P- Moreover, the optimal control is constant and given 
by : 



aw 
~dt 

dw 



[t,x) 
(t,x) 
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The key point in the explicit resolution of the Merton problem is that the value 
function v may be separated into a function of t and of x :v(t,x) = tp(t)x p . With 
this transformation and substitution into the HJB equation, it turns out that 
if is the solution of an ordinary differential equation with terminal condition 
ip(T) = 1, which is explicitly solved. Other applications of verification theorems 
to stochastic control problems in finance are developed in the recent textbook 
|7U| . There are also important examples of applications in engineering with 
the stochastic linear regulator, which will be studied later by the maximum 
principle. 

Classical verification theorems allow to solve examples of control problems 
where one can find, or at least there exists, a smooth solution to the associated 
HJB equation. They apply successfully for control problems where the diffusion 
term does not depend on the control and is uniformly elliptic, since in this 
case the HJB equation is semilinear in the sense that D\v appears linearly, 
and so classical existence results for smooth solutions exist, see [HO]. They also 
apply for some specific models with control on the diffusion term, typically 
Merton's portfolio selection problem as described above, and more generally to 
extensions of Merton's model with stochastic volatility, since in this case, the 
HJB equation may be reduced after a suitable transformation to a semilinear 
equation. This last point is developed in [51] and |77]. However, in the general 
case of nondegeneracy of the diffusion term and in a number of applications, the 
value function might be not smooth or it is not possible to obtain a priori the 
required regularity. Moreover, for singular control problems, the value function 
is in general not continuous at the terminal date T, so that the right terminal 
condition is not given by g(x), i.e. lim f /Tii(f,i) ^ g(x). Then, the classical 
verification approach does not work and we need to relax the notion of solution 
to the HJB equation. It turns out that the suitable class of weak solutions is 
the one of viscosity solutions, not only for dealing with the rigorous derivation 
of the HJB equation, but also for determining the relevant terminal condition. 
Moreover, viscosity solutions theory provides a general verification approach, 
which allows to go beyond the case of classical verification. We shall give some 
applications in finance of the viscosity solutions approach where we can explicit 
non smooth solutions to singular control problems. 

The notion of viscosity solutions was introduced by P.L. Lions [HI] for second 
order equations. They provide by now a well-established method for dealing with 
stochastic control problems, in particular in finance, see e.g. [50], [H], (EPJ or 
|5T] . We review some of these aspects and their applications in the next section. 

2.3. Pontryagin's maximum principle 

A classical approach for optimization and control problems is to derive necessary 
conditions satisfied by an optimal solution. For example, the Kuhn- Tucker condi- 
tion is a well-known necessary condition of optimality in the finite-dimensional 
case. For an optimal control problem as described above, and which may be 
viewed as an infinite-dimensional optimization problem, the argument is to use 
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an appropriate calculus of variations on the gain function J(t, x, a) with respect 
to the control variable in order to derive a necessary condition of optimality. 
The maximum principle, initiated by Pontryagin in the 1960s, states that an 
optimal state trajectory must solve a Hamilton system together with a maxi- 
mum condition of a function called a generalized Hamilton. In principle, solve 
a Hamilton should be easier than solving the original control problem. 

The original version of Pontryagin's maximum principle was derived for de- 
terministic problems. As in classical calculus of variation, the basic idea of is to 
perturb an optimal control and to use some sort of Taylor expansion of the state 
trajectory and objective functional around the optimal control. By sending the 
perturbation to zero, one obtains some inequality, and by duality, the maximum 
principle is expressed in terms of an adjoint variable (Lagrange multiplier in the 
finite-dimensional case). The stochastic control case was extensively studied in 
the 1970s by Bismut, Kushner, Bensoussan or Haussmann. However, at that 
time, the results were essentially obtained under the condition that there is no 
control on the diffusion coefficient. For example, Haussmann investigated max- 
imum principle by Girsanov's transformation and this limitation explains why 
this approach does not work with control-dependent and degenerate diffusion 
coefficients. The main difficulty when facing a general controlled diffusion is that 
the Ito integral term is not of the same order as the Lebesgue term and thus the 
first-order variation method fails. This difficulty was overcomed by Peng |75| . 
who studied the second-order term in the Taylor expansion of the perturbation 
method arising from the Ito integral. He then obtained a maximum principle for 
possibly degenerate and control-dependent diffusion, which involves in addition 
to the first-order adjoint variable, a second-order adjoint variable. 

In order to make applicable the maximum principle, one needs some explicit 
description of the adjoint variables. These variables obtained originally by dual- 
ity in functional analysis may be represented by Riesz representation of a certain 
functional. By completing with martingale representation in stochastic analy- 
sis, the adjoint variables are then described by what is called today backward 
stochastic differential equations (BSDE). Actually, it was the study of the max- 
imum principle in the stochastic control case that motivated Peng for a general 
formulation of BSDE, which in turn generated an important area of research. 
We shall state in Section 0] the mathematical formulation of the sufficiency of 
the maximum principle and relate it to BSDE and its extensions. 

2-4- Other control problems 

We present in this paragraph some other control problems, which we do not 
study in detail here, but that are also of significant theoretical and practical 
interest. We also emphasize some present developments. 

2. If.. 1. Random horizon 

In problem formulation (|2.2|) . the time horizon is fixed until a deterministic 
terminal time T. In some real applications, the time horizon may be random, 
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and in the context of controlled diffusion described in paragraph l2.il the control 
problem is formulated as : 



where r is a finite random time. In standard cases, the terminal time r is a 
stopping time at which the state process exits from a certain relevant domain. 
For example, in a reinsurance model, the state process X is the reserve of a 
company that may control it by reinsuring a proportion 1 — a of premiums 
to another company. The terminal time r is then the bankruptcy time of the 
company defined as r = inf{i > : X t < 0}. More generally, given some open 
set O of R n , 



(which depends on the control). In this case, the control problem H2.11(l leads 
via the dynamic programming approach to a Dirichlet boundary- value problem. 
Another case of interest concerns a terminal time r, which is a random time but 
not a stopping time in the filtration F with respect to which the controls are 
adapted. This situation occurs for example in credit risk models where r is the 
default time of a firm. Under the so-called (G) hypothesis on filtration theory : 
P[t < t\T t ] is a nondecreasing right-continuous process, problem p.llj) may 
be reduced to a stochastic control problem under a fixed deterministic horizon, 
see |13 | for a recent application in portfolio optimization model. In the general 
random time case, the associated control problem has been relatively lightly 
studied in the literature, see |17| or |92| for a utility maximization problem in 
finance. 

2-4-2. Optimal stopping 

In the models presented above, the horizon of the problem is either fixed or indi- 
rectly influenced by the control. When one has the possibility to control directly 
the terminal time, which is then modelled by a controlled stopping time, the 
associated problem is an optimal stopping time problem. In the general formu- 
lation of such models, the control is mixed, composed by a pair control/stopping 
time (a, r) and the functional to optimize is : 



The theory of optimal stopping, thoroughly studied in the seventies, has received 
a renewed interest with a variety of applications in economics and finance. These 
applications range from asset pricing (American options) to firm investment 
and real options. Extensions of classical optimal stopping problems deal with 
multiple optimal stopping with eventual changes of regimes in the state process. 
They were studied e.g. in [5], [87], and applied in finance in PJ], [301, 05], |2T| 
or ISSj). 




(2.11) 



r 



inf {t > : X t i O} A T 
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2.4-3. Impulse control 

In formulation of the problem in paragraph 12.11 the displacement of the state 
changes continuously in time in response to the control effort. However, in many 
real applications, this displacement may be discontinuous. For example, in in- 
surance company models, the company distributes the dividends once or twice 
a year rather than continuously. In transaction costs models, the agent should 
not invest continuously in the stock due to the costs but only at discrete times. 
A similar situation occurs in a liquidity risk model, see e.g. |2~2" . Impulse con- 
trol provides a suitable framework for modelling such situations. This may be 
described as follows : the controlled state diffusion process is governed by 



where the control C is a pure jump process. In other words, the control is given 
by a pair (t„,k„)„ where (r„)„ is a nondecreasing sequence of stopping times, 
representing the intervention times of the controller, and {n n ) n is a sequence 
of T Tn -measurable random variables, representing the jump size decided by the 
controller at time r„. The functional objective to optimize is in the form : 



Impulse control problem is known to be associated via the dynamic program- 
ming approach to an HJB quasi- variational inequality, see For some recent 
applications in finance, we refer to for insurance models, [HHj and pHl fc> r 
transaction costs models, and for liquidity risk model. 

2-4-4- Ergodic control 

Some stochastic systems may exhibit over a long period a stationary behavior 
characterized by an invariant measure. This measure, if it does exists, is obtained 
by the average of the states over a long time. An ergodic control problem consists 
in optimizing over the long term some criterion taking into account this invariant 
measure. 

A standard formulation resulting from the criterion presented in paragraph 
12. H is to optimize over control a functional of the form : 



dX s 



b{s, X s )dt + ct(s, X s )dW s + dC 





or 
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This last formulation is called risk-sensitive control on an infinite horizon. Er- 
godic and risk-sensitive control problems were studied in JT] or Risk- 
sensitive control problems were recently applied in a financial context in |12| 
and [%U| . 

Another criterion is based on the large deviations behavior of the ergodic 
system : P[Xt/T] ~ e _/( - c ^ T , when T goes to infinity, consists in maximizing 
over control a functional of the form : 

lim sup — In P 

This large deviations control problem is interpreted in finance as the asymptotic 
version of the quantile criterion of maximizing the probability that the terminal 
wealth Xt beats a given benchmark. This nonstandard control problem was 
introduced and developed recently by Pham [7S| , [75| , see also @7] . It does not 
have a direct dynamic programming principle but may be reduced via a duality 
principle to a risk-sensitive control problem. 

2.4-5. Robust control 

In the problems formulated above, the dynamics of the control system is assumed 
to be known and fixed. Robust control theory is a method to measure the 
performance changes of a control system with changing system parameters. This 
is of course important in engineering systems, and it has recently been used in 
finance in relation with the theory of risk measures initiated by Artzner et al. 

Indeed, it is proved that a coherent risk measure for an uncertain payoff Xt 
at time T is represented by : 

p(-X T ) = sup E Q [X T ], 
QeQ 

where Q is a set of absolutly continuous probability measures with respect to 
the original probability P. More generally, one may define a risk measure in the 
form : 

p(-X T ) = - mtEQ[U(-X T )\, 

where U is a concave and nondecreasing function. So, when X is controlled by a 
and the problem is to optimize the risk measure p(—Xx), one is facing a robust 
control problem. In this financial context, robust optimization problems were 
recently studied in (221 and 

2-4-6. Partial observation control problem 

It is assumed so far that the controller completely observes the state system. In 
many real applications, he is only able to observe partially the state via other 



Xt 



> c 
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variables and there is noise in the observation system. For example in financial 
models, one may observe the asset price but not completely its rate of return 
and/or its volatility, and the portfolio investment is based only on the asset 
price information. We are facing a partial observation control problem. This 
may be formulated in a general form as follows : we have a controlled signal 
(unobserved) process governed by 

dX s = b(s,X s ,Y s ,a s )ds + a(s,X s ,Y s ,a s )dW s , 

and an observation process 

dY s = r](s,X s ,Y s ,a s )ds + i(s,X s ,Y s ,a s )dB s , 

where B is another Brownian motion, eventually correlated with W. The control 
a is adapted with respect to the filtration generated by the observation ¥ Y = 
(Jt) an d the functional to optimize is : 



J(a) = E 



I f(X t ,Y u a t )dt + g(X T ,Y T ) 
Jo 



By introducing the filter measure-valued process 

n t (da;) = P[X t €dx\7?], 
one may rewrite the functional J (ex) in the form : 



J(a) 



E 



f(H t ,Y t ,a t )dt + g(IL T ,Y T ) 



where we use the notation : f(ir,y) = J f(x,y)i:(dx) for any finite measure 7r 
on the signal state space, and similarly for g. Since by definition, the process 
(n t ) is (J r t y )-adapted, the original partial observation control problem is refor- 
mulated as a complete observation control model, with the new observable state 
variable defined by the filter process. The additional main difficulty is that the 
filter process is valued in the infinite-dimensional space of probability measures : 
it satisfies the Zakai stochastic partial differential equation. The dynamic pro- 
gramming principle or maximum principle are still applicable and the associated 
Bellman equation or Hamiltonian system are now in infinite dimension. For a 
theoretical study of optimal control under partial observation under this infinite 
dimensional viewpoint, we mention among others the works |36) . |26j . |l(Jj . 
1^1] or There are relatively few explicit calculations in the applications to 
finance of partial observation control models and this area should be developed 
in the future. 



2.4-7. Stochastic target problems 



Motivated by the superreplication problem in finance, and in particular under 
gamma constraints |HS], Soner and Touzi introduced a new class of stochastic 
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control problems. The state process is described by a pair (X, Y) valued in 
R™ x R, and controlled by a control process a according to : 

dX s = b(s,X s ,a s )ds + a(s,X s ,a s )dW s (2.12) 
dY s = i7(a,jr a ,y;,a a )d8 + 7(s,JC s ,y s ,a a )dW.. (2.13) 

Notice that the coefficients of X do not depend on Y. Given a E A and (t, x, y) 
S [0, T] x R" x R, (X*'*, Y*'**) is the unique solution to (f?TT^) - (f^TT5f> with initial 
condition (X t ' x , Y^' x,v ) = (x, y). The coefficients 6, <r, rj, 7 are bounded functions 
and satisfy usual conditions ensuring that (X f - X , Y t ' x ' v ) is well-defined. The 
stochastic target problem is defined as follows. Given a real- valued measurable 
function g on R™ , the value function of the control problem is defined by : 

v(t,x) = infjyeR :3a e A, Y^' x,y > g(X^ x ) a.s.} . 

In finance, X is the price process, Y is the wealth process controlled by the 
portfolio strategy a, and v(t, x) is the minimum capital in order to superreplicate 
the payoff option q{Xt)- 

The dynamic programming principle associated to this stochastic target prob- 
lem is stated as follows : for all (t, x) 6 [0, T], and 9 stopping times in [t, T], we 
have 

v(t,x) = infjyeR :3aeA, Y*< x ' y > u(6>,X*' x ) a.s.}. 

The derivation of the associated dynamic programming equation is obtained 
under the following conditions. The matrix a(t, x, a) is assumed to be invertiblc 
and the function 

a 1 > jit, x, y, a)cr _1 (i, x, a) 

is one-to-one for all (t,x,y), with inverse denoted 1?, i.e. 

7(t,x, y,a)cr _1 (t,a;, a) = p' a = i3(t,x,y,p) 

for p E R". Here p' is the transpose of p. Moreover, the control set A is assumed 
to be convex and compact in R m , with nonempty interior. The support function 
of the closed convex set A is denoted 

MC) = supa'C- 

aeA 

Under these conditions, Soner and Touzi |8fi| proved that v is a viscosity solution 
to the dynamic programming equation : 

min x) - C ao(t ' x) v(t, x) + T)(t, x, v(t, x), a (t, x)) , 

G(t,x,v(t,x),D x v(t,x))} = 0, 

where C a is the second order differential operator associated to the diffusion X, 
and 

ao(t,x) = tf(t,x,v(t,x),D x v(t,x)), 
G(t,x,v(t,x),D x v(t,x)) = inf [5 A (() ~ a (t,x)'Q ■ 
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3. Dynamic programming and viscosity solutions 
3.1. Definition of viscosity solutions 

The notion of viscosity solutions provides a powerful tool for handling stochastic 
control problems. The theory of viscosity solutions goes beyond the context of 
HJB equations and is a concept of general weak solutions for partial differential 
equations. We refer to the user's guide of Crandall, Ishii and Lions for an 
overview of this theory. Here, we simply recall the definition and some basic 
properties required for our purpose. 

Let us consider parabolic nonlinear partial differential equations of second 
order : 

F(t,x,w(t,x),^-(t,x),D x w(x),D xx w{x)) = 0, (t, x) e [0, T) x 0(3.1) 
at 

where is an open set of R™, and F is a continuous function on [0, T] x x R x 
R x R™ x S n . The function F is assumed to satisfy the ellipticity and parabolicity 
conditions : 

M<M F(t,x,r,p t ,p,M) > F(t,x,r, Pt ,p,M) (3.2) 

Pt<Pt => F(t,x,r, Pt ,p,M) > F(t,x,r,p t ,p,M), (3.3) 

for all t e [0,T), x G O n , r £ R, p t ,p t 6 1, p e 1 and M,M G «S„. The last 
condition l|3.3l) means that we are dealing with forward PDE, i.e. (13.11) holds for 
time t < T, and the terminal condition is for t = T . This is in accordance with 
the control problem and HJB equation formulated in Section [3 

Since it is not always easy to have a priori the continuity of the value function 
(actually, there may be even discontinuity at terminal time T), we work with 
the notion of discontinuous viscosity solutions. We then introduce, for a locally 
bounded function w on [0, T] x 0, its lower semicontinuous envelope w* and 
upper semicontinuous envelope w* , i.e. : 

w*(t, x) = liminf w(t',x') and w*(x) = limsup w(t',x'). 

(*'.»')-•<*»«) [«',»')-.(t,») 

Definition 3.1 Let w be a locally bounded function on O. 

(i) w is a viscosity subsolution (resp. supersolution) of (|3.1(l if : 

F{t,x,ip{i,x)^(i,x),D xl p{t,x),Dl<p{t,x)) < (resp. >) 0, (3.4) 

for any (t, x) G [0, T) x and smooth test function ip G C 2 ([0, T) x 0) s.t. (t, x) 
is a maximum (resp. minimum) ofw* — (p (resp. W* — (p) withO = (w* — tp)(t,x) 
(resp. (to, - <p)(i,x)). 

(ii) w is a viscosity solution of ()3.1|l if it is a viscosity subsolution and superso- 
lution. 
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Remark 3.1 1. Viscosity solutions extend the notion of classical solutions : a 
C 1,2 function on [0, T) x O is a supersolution (resp. a subsolution) of (|3.1I) in 
the classical sense iff it is a viscosity supersolution (resp. viscosity subsolution). 

2. The above definition is unchanged if the maximum/minimum is strict and/or 
local. Notice that since we are considering forward PDEs, a local extremum 
point (t, x) means with respect to a neighborhood of the form [t, t + h) x B n (x), 
with h,rj > 0. Here B n (x) denotes the open ball of radius r\ and center x and 
B v {x) its closure. 

3. Without loss of generality, by translating the test function by a constant, wc 
may relax the condition that = (w* — (p)(i,x) (resp. (w* — cp)(i,x)). We have 
then to replace in (|3.4(l ip(i, x) by w*(i,x) (resp. w*(i, x)). 

4. We define similarly viscosity solutions for elliptic PDEs, i.e. without time 
variable t. 

3.2. Viscosity characterization 

We come back to the framework of controlled diffusions of Section^ and we state 
that the value function is a viscosity solution to the associated HJB equation. 
We also determine the relevant terminal condition. We present a unifying result 
for taking into account the possible singularity of the Hamiltonian H when the 
control set A is unbounded, see Remark 12.21 We then introduce 

dom(H) = {(t, x,p, M) e [0, T] x E" x M™ x S n : H(t, x,p, M) < 00} , 

and we shall make the assumption : 

H is continuous on int(dom(i/)) 
and there exists G : [0, T] x R" x M™ x S n 
nondecreasing in its last argument and continuous s.t. : 

(t,x,p,M) edom(H) <=^> G(t,x,p,M)>0 (3.5) 

The nondecreasing condition of G in its last argument is only required here to 
ensure that the PDE in (|3.6|) satisfies the ellipticity condition (|3.2|l . 

In the sequel, we assume that the value function v is locally bounded on 
[0, T] x R" : this is a condition usually easy to check, and satisfied typically 
when v inherits from / and g a linear growth condition. Our general viscosity 
property for the value function is stated in the following theorem : 

Theorem 3.1 Assume that /(., ., a) is continuous for all a € A and Ij3.5|l holds. 
Then, v is a viscosity solution of the HJB variational inequality : 

{dv 
~ ~di ^' ^ _ H ^' X ' DxV ( t ' D ^ v ^ 7 x ^ ' 

G{t, x, D x v{t, x), Dlv(t, x))} = 0, (t,x)E [0,T) x R". (3.6) 
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Remark 3.2 In the regular case, i.e. when the Hamiltonian H is finite on the 
whole state domain (this occurs typically when the control set A is compact), 
the condition l|3.5|) is satisfied for any choice of positive continuous function, e.g. 
a positive constant. In this case, the HJB variational inequality l|3.6[) is reduced 
to the regular HJB equation : 

-—(t,x)-H(t,x,D x v(t,x),Dlv(t,x)) = 0, ((,i)e[0,r)xl" 
at 

that the value function satisfies in the viscosity sense. Hence, Theorem l3.1l states 
a general viscosity property including both the regular and singular case. 

Remark 3.3 We do not address here the important uniqueness problem asso- 
ciated to the HJB (in)equation. We refer to [231 f° r some general uniqueness 
results. In most cases, there is a (strong) comparison principle for this HJB 
PDE, which states that any uppersemicontinuous subsolution is not greater 
than a lowersemicontinuous supersolution. This implies v* < u* and therefore 
v* = since the other inequality is always true by definition. The consequence 
of this is the continuity of the value function v on [0, T) x K n . Hence, we notice 
that the viscosity property allows to derive as a byproduct the continuity of the 
value function, which may not always be easily proved by a direct probabilistic 
argument. 

It is well known that to a parabolic PDE is associated a terminal condition, 
in order to get a uniqueness result. We then need to determine the terminal data 
for the value function. By the very definition of the value function, we have 

v(T,x) = g(x), xe~R d . (3.7) 

However, in several aplications of stochastic control problems, the value function 
may be discontinuous at T, see e.g. [20] and |23]. In this case, (|3.7(l is not 
the relevant terminal condition associated to the HJB equation in order to 
characterize the value function : we need actually to determine v(T~,x) :— 
\im t yTv(t,x) if it exists. 
To this end, we introduce 

v(x) = liminf v(t,x') and v(x) — limsup v(t,x') 

t/T,x'^,x tST,x'->x 

Notice that by definition, v < v, v is lowersemicontinuous, and v is uppersemi- 
continuous. 

Our characterization result for the terminal condition is stated as follows. 

Theorem 3.2 Assume that f and g are lower-bounded or satisfy a linear growth 
condition, and (|3.5|l holds. 

(i) Suppose that g is lowersemicontinuous. Then v is a viscosity supersolution 
of 

min \_v_{x ) -g(x) , G{T,x,Dv{x),D 2 v(x))] = 0, on 1". (3.8) 
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(ii) Suppose that g is uppersemicontinuous. Then v is a viscosity subsolution of 

min [v{x) -g(x) , G{T,x,Dv{x) 7 D 2 v{x))] = 0, on R n . (3.9) 

Remark 3.4 In most cases, there is a comparison principle for the PDE arising 
in the above theorem, meaning that a subsolution is not greater than a superso- 
lution. Therefore, under the conditions of Theorem 13. 21 we have v < v and so v 
= v. This means v := v(T~ , .) exists, equal to v — v and is a viscosity solution 
to : 

min \v(x) — g(x) , G(T,x,Dv(x),D 2 v(x))] = 0, on R n , (3.10) 

Denote, by g the upper G-envelope of g, defined as the smallest function above 
g and viscosity supersolution to : 

G{T,x,Dg{x),D 2 g{x)) = 0, on M™, (3.11) 

when it exists and is finite. Such a function may be calculated in a number of 
examples, see e.g. paragraph s. 3. 21 Since v is a viscosity supersolution to (|3.1U|) . 
it is greater than g and is a viscosity supersolution to the same PDE as g. Hence, 
by definition of g, we have v > g. On the other hand, g is a viscosity supersolution 
to the PDE H3.1()|l . and so by a comparison principle, the subsolution ii of H3.10[l 
is not greater than g. We have then determined explicitly the terminal data : 

v(T~,x) = v(x) = g(x). 

Recall that in the regular case, we may take for G a positive constant function, 
so that obviously g = g. Therefore, in this case, v is continuous in T and w(T _ , x) 
= v(T, x) = g(x). In the singular case, g is in general different from g and so v is 
discontinuous in T. The effect of the singularity is to lift up, via the G operator, 
the terminal function g to g. 

We separate the proof of viscosity supersolution and subsolution, which are 
quite different. The supersolution part follows from DPI and standard argu- 
ments in viscosity solution theory as in Lions 61 . In this part, the notion of 
viscosity solutions is only used to handle the lack of a priori regularity of v. 
The subsolution part is more delicate and should take into account the possible 
singular part of the Hamiltonian. The derivation is obtained from DP2 and a 
contraposition argument. The assertion of the subsolution part as well as the 
determination of the terminal condition seem to be new in this unifying context 
and are inspired by arguments in |8f>j . These arguments really use the concept 
of viscosity solutions even if the value function were known to be smooth. The 
reader who is not interested in the technical proofs of Theorems 13.11 and 13.21 
can go directly to paragraph 13.31 for some applications in finance of viscosity 
solutions theory. 
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3.2.1. Viscosity supersolution property 

Proposition 3.1 Assume that /(., ., a) is continuous for all a £ A. Then, v is 
a viscosity supersolution of the HJB equation : 

-—(t,x)-H(t,x,D x v(t,x),Dlv(t,x)) = 0, (t, x) £ [0, T) x R™. (3.12) 

Proof. Let (i, x) £ [0, T) x R" and tp £ C 2 ([0, T) x R n ) a smooth test function 
s.t. : 



= («* - ip)(t,x) 



min (u* — <p)(t, x). 
((,i)e[o,r)xr 



(3.13) 



By definition of v*(t,x), there exists a sequence (t m ,x m ) in [0, T) x R" s.t. 

(*m,aj m ) {i,x) and u(t m ,ar m ) v*(t,ac), 
when m goes to infinity. By continuity of <p and H3.13[) . we also have 

when to goes to infinity. 

Let a be an arbitrary element in A, and a the constant control equal to a. 
We denote by X t s m ' Xm the associated controlled process starting from x m at t m . 
Consider r m the first exit time of X tm ' Xm from the open ball B„(x m ) : r m = 
inf{s > t m : \Xl m ' Xm — x m \ > 77}, with 77 > and let (h m ) a positive sequence 
s.t. : 



and 

' ' /// 



0. 



when to goes to infinity. By applying the first part (DPI) of the dynamic pro- 
gramming principle to v(t m ,x m ) and m := r m A (t m + h m ), we get : 



v(t m ,x m ) > E 



f(s, xy xm ,a)ds + v(e m ,xl2' Xm ) 



From l|3.13|l . which implies v > v* > tp, we then deduce : 



) + j m > E 



i)ds + <p(9 m , xl m ' Xm ) 



We now apply Ito's formula to <p(s, Xl m ' Xm ) between t m and 9 m , and we obtain 
after noting that the stochastic integral term vanishes in expectation due to 
bounded integrand : 



7m 
hrn 



E 



dtp 
~dt 



C a p-f)(s,X t - x -,a)ds 



> 0. (3.14) 
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By the a.s. continuity of the trajectory X l s m,Xm , we get that for m sufficiently 
large, (m > N {(*))), 9 m {uo) — t m + h m , a.s. We then deduce by the mean- value 
theorem that the random variable inside the expectation in (|3.14|) converges a.s. 
dtp 

to — — (t,x) — C a (p(i,x) — f(t,x,a) when m goes to infinity. Moreover, since 

this random variable is bounded by a constant independent of m, we may apply 
the dominated convergence theorem to obtain : 

-^(i,x)-£°<p(i,x)-f(t,x,a) > 0. 

We conclude from the arbitrariness of a G A. □ 

Remark 3.5 The supersolution property of v means that for all (t, x) <E [0, T) x 
R™ and smooth test function (p s.t. = (u* — tp)(i, x) = min[o ! T)xR"( t; * — <p): we 
have 

Qu _ 

-—(i,x)-H{i,x,D x (p(i,x),Dl<p(t,x)) > 0. 

Recalling condition (|3.5() . this implies 

G(i,x,D x ip(i,x),D 2 x ip(i,x)) > 0, 

and so 

~ — {t,x) - H(i,x,D x tp{i,x),Dl(p(t,x)) , 
G(t,x,D x <p(i,x),Dltp(i,x))} > 0. 
This is the viscosity supersolution of v to (|3.6() . 



3.2.2. Viscosity subsolution property 

Proposition 3.2 Assume that Ij3.5|l holds. Then v is a viscosity subsolution of 
the HJB variational inequality : 

{du 
- — (t, x) ~ H(t, x, D x v(t, x),D x v(t, x)) , 

G{t, x, D x v{t, x), D 2 x v{t, x))} = 0, (t,x)E [0,T) x R".(3.15) 

The proof of the subsolution part is based on a contraposition argument and 
DP2. We then introduce for a given smooth function ip, the set in [0, T] x R n : 

M(tp) = { (t, x) e [0, T] x R™ : G{t, x, D x ip(t, x), D x ip(t, x)) > 

and ~ -^-(t,x) ~ H(t,x,D x ip(t,x),D 2 x ip(t,x)) > 

The following Lemma, which will be also used in the derivation of the terminal 
condition, is a consequence of DP2 of the dynamic programming principle. 
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Lemma 3.1 Let ip be a smooth function on [0, T] x R™, and suppose there exist 
h < t 2 < T, x G R n and r, > s.t. : 

[ti,t 2 ] x B v (x) e M(<p). 

Then, 

sup (v — tp) = max (v* — tp), 

d p ([t 1 ,t 2 ]xB v (x)) [ti,t 2 ]xB n (S) 

where d p {[ti,t2\ x B v (x)) is the forward parabolic boundary of [tijta] x B n {x), 
i.e. d p ([ti,t 2 ] x B„{x)) = [ti,t 2 ] x SB, (as) U {< 2 } x B v (x). 

Before proving this Lemma, let us show how it immediately implies the re- 
quired subsolution property. 

Proof of Proposition l3.21 Let (t, x) £ [0, T) xl" and tp a smooth test function 
s.t. 

= (v* — tp)(i,x) = (strict) max (v* — tp). 

[0,T)xR" 

First, observe that by the continuity condition in Q3.5|l . the set M.(tp) is open. 
Since (t, x) is a strict maximizer of (v* — tp), we then deduce by Lemma l3.1l that 
(i,x) Ai(ip). By definition of M(tp), this means : 



mm 



~{t,x) -H(t,x,DMt,x),Dl<p(t,x)) , 



G(t,x,D x tp{t,x),D 2 x tp{t,x))} < 0, 
which is the required subsolution inequality. 

Proof of Lemma 13.11 By definition of A4(tp) and H, we have for all a e A : 

- ^(t,x) -C a tp{t,x) -f(t,x,a) > 0, V(t,x) G [t x ,t 2 ] x 5„(i)(3.16) 
We argue by contradiction and suppose on the contrary that : 

max (v* — tp) — sup (v — tp) := 28. 

[t 1 ,t 2 ]xB, 1 (x) ' d p ([t 1 ,t 2 ]xB v (x)) 

Wecan choose (to,xo) G {t\,t 2 )Y.B r ,(x) s.t. (v—tp)(to,xo) > —6+ max (v*- 

[t 1 ,t 2 ]xB v (x) 

tp), and so : 

(v - tp)(t ,x ) > 6+ sup (v-tp). (3-17) 

d p {[t 1 ,t 2 \xB rl (x)) 



Fix now e = 6/2, and apply DP2 to v(to, xq) : there exists a £ £ A s.t. 



u(to) ^o) — £ < -E 



(3.18) 



H. Pham/On some recent aspects of stochastic control and their applications 528 



where we choose 



inf {s > t Q : (s, X^ x °) $ [h,t 2 ] x B v (x)} . 



First, notice that by continuity of X to > x ° , we have (9,X g °' x °) G d p {[t 1 ,t 2 ] x 
B n (x)). Since from l|3.17[) . we have v < <p + (v — ip)(t ,xo) — S on d p ([ti,t 2 ] x 
Br,(x)), we get with lj3~T%|) : 



f(s, Xl°' x °,a e s )ds + p(9, Xl°' X0 ) - tpfa, so) 



5. 



Applying Ito's formula to y>(s, JT* ' 10 ) between s = t and s = 9, we obtain : 



in 



e-S. 



Since, by definition of 9, (s,Xl°> x °) lies in [ti,t 2 ] x B v (x) for all t < s < 9, we 
get with (|3.16(l the required contradiction : < e — S — —5/2. 

3.2.3. Terminal condition 



We start with the following Lemma. 

Lemma 3.2 Suppose that f and g are lower-bounded or satisfy a linear growth 
condition, and g is lower semicontinuous. Then, 

v{x) > g(x), Vx G 1". 

Proof. Take some arbitrary sequence (t mi x m ) — > (T,x) with t m < T and fix 
some constant control a = a G A. By definition of the value function, we have : 



> E 



f(s 1 X t s - x -,a)d S + g(X t T - x -) 



Under the linear growth or lower-boundeness condition on / and <?, we may 
apply the dominated convergence theorem or Fatou's lemma, and so : 

liminf v(t m , x m ) > E liminf g(X^ l ' Xrr 

by the lowersemicontinuity of g and the continuity of the flow Xl^ x in (t, x). □ 

The supersolution property l|3.8|) for the terminal condition is then obtained 
with the following result. 

Lemma 3.3 Under J33j, v is a viscosity supersolution of : 
G(T,x,D x v(x),D 2 x v{x)) = 0, on M™. 
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Proof. Let x £ W 1 and ip a smooth function on R" s.t. 

= (v-ip)(x) = min(v - ip). (3.19) 

By definition of v, there exists a sequence (s m , y m ) converging to (T, x) with s m 
< T and 

lim v*(s m ,z/ m ) = v(x). (3.20) 

m — >oc 

Consider the auxiliary test function : 

T — t 

<Pm(t,x) = ljj(x) - \X - x\ 4 + — rj, 

and choose (t m , z m ) £ [s m , T] x B\(x) as a minimum of (v* — ^? m ) on [s m , T] x 
Bi(x). 

5£ep 1. We claim that, for sufficiently large m, t m < T and x m converges to 
x, so that (t m ,x m ) is a local minimizer of (i>* — <p m )- Indeed, recalling vjx) — 
il>(x) and Q3.2Up . we have for sufficiently large m : 

(«* - <p m )(s m ,y m ) < -rr^r- V < °- ( 3 - 21 ) 

2(T - s m ) 

On the other hand, for any x £ K n , we have 

(«* -y» m )(T,x) = v{x) - i)(x) + \x - x\ A > v(x)-ip(x) > 0,(3.22) 

by (jXTOjl . The two inequalities (|3~2T|) - (|3~22|) show that t m < T for large m. We 
can suppose that x m converges, up to a subsequence, to some Xq £ B\(x). From 
<|3.19[1 , since s m < t m and (i m , x m ) is a minimum of (v* — ip m ) , we have : 

< (v-i/j)(x )-(v-iP)(x) 

< liminf Uv* - tp m )(t m ,x m ) - (u* - ip m )(s m ,y m ) - \x m - x\ A ] 

m—*oo 

< -Ixo-xl 4 , 

which proves that xq = x. 

Step 2. Since (t m ,x m ) is a local minimizer of (v* — (p m ), the viscosity super- 
solution property of v* holds at (t m ,x m ) with the test function tp m , and so for 
every m : 

G (tm , X m , D x (p m (tm i ) : D x (p m (t m : %m ) ) — 0. (3.23) 

Now, since D x (p m (t m ,x m ) = Dip(x m ) - A(x m - x)\x m - x\ 2 , D 2 (p m (t m ,x m ) = 
D 2 ip(x m ) — 4|x TO — x\ 2 I n — 4(x m — x)(x m — x)' , recalling that G is continuous, 
and (t m ,x m ) converges to (T,x), we get from (|3.23l) : 

G(T,x,DiP(x),D 2 i>{x)) > 0. 

This is the required supersolution inequality. □ 

We finally turn to the subsolution property for the terminal condition. 



H. Pham/On some recent aspects of stochastic control and their applications 530 

Lemma 3.4 Suppose that g is uppers emicontinuous and (|3.5|1 holds. Then, v 
is a viscosity subsolution of : 

min [v(x) — g(x) , G(T, x, Dv(x), D 2 v(x))] = 0, on K" 

Proof. Let x £ K™ and tp a smooth function on l n s.t. 

= (v-tp)(x) = max(v-ip). (3.24) 

We have to show that whenever 

v(x) > g(x), (3.25) 

then 

G{T,x,DiP(x),D 2 iP(x)) < 0. (3.26) 
So, suppose that H3.25JI holds, and let us consider the auxiliary test function : 

(f m (t,x) = ip(x) + \x — x\ A + m(T — t). 
We argue by contradiction and suppose on the contrary that 

G(T,x,DiP(x),D 2 ij(x)) > 0. 

Since D x ip m (t, x) = Dip(x) —A(x — x)\x — x\ 2 — * Difj(x), D 2 ip m (t, x) = D 2 tp(x) — 
4:I n \x — x\ 2 — 4(x — x){x — x)' —¥ D 2 ip(x) when x tends to x, there exists, by 
continuity of G, so < T and r\ > s.t. for all m : 

G(t, x, D xVm {t, x), D 2 x p m (t, x)) > 0, V(f, x) £ [s , T] x B n (x). (3.27) 

Under condition H3.5(l . the function H (t, x, D x ip m (t, x), D 2 (p m (t, x)) is then finite 
on the compact set [so, T] x B n (x) and by continuity of H on int(dom(i7)), there 
exists some constant ho (independent of m) s.t. 

H(t,x,D x ip m (t,x),Dlip m (t,x)) < h Q , V(t,x) £ [so,T] x B v {x). (3.28) 

Step 1. Since by definition v > v, we have from Lemma 13.21 : 

v > g. (3.29) 

Hence, for all x £ E n , 

(v - tp m )(T, x) = (g - ip)(x) - \x - x\ A < (v - tp)(x) - \x - x\ 4 

< ~\x-x\ 4 < (3.30) 

by (|3.24|) . This implies : swp Bfi ^(v — ip m )(T,.) < 0. We claim that 

limsup sup (v — ip m )(T, .) < 0. (3.31) 

■m^oo B rl (x) 
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On the contrary, there exists a subsequence of (<p m ), still denoted ((p m ) s -t- : 
lim sup (v — <p m )(T, .) = 0, 

For each m, let be a maximizing sequence of (v — ip m )(T, .) on B v (x), i.e. 

lim lim («-¥> m )(T,a4) = 0. 

Now, from l|3.3U|) . we have (w — (f m )(T, x^) < — \x m ~ ^| 4 > which combined with 
the above equality shows that : 



lim lim x m = x. 

m^oo k— »oc 



Hence, 



= lim lim (v — ip m )(T, ajjjj = lim lim g{x k m ) — ip(x) 
< g(x) - i/j(x) < (v - %/})(x), 
by the uppersemicontinuty of g and (|3.25|) . This contradicts (v — ip)(x) = in 

iJXHJi. 

Step 2. Take a sequence (s m ) converging to T with s < s m < T. Let us 
consider a maximizing sequence (i m , x m ) of v* — ip m on [s m , T] x dB v (x). Then 

limsup sup (w* - <yj m ) < limsup(w*(t m , x m ) - ip(x m )) ~ V 4 - 

m ~>°° [s m ,T]xdB v (x) m-^oo 

Since t m converges to T and x m , up to a subsequence, converges to some xo € 
dB v {x), we have by definition of "0 : 

limsup sup (y* — y> m ) < (v — i())(xo) — r/ 4 < — r? 4 , (3.32) 

m ^°° [s m ,T]xdB v (x) 

by l|3~2^jl . Recall also from (j3~2I|) that (v* - <p m )(T, x) = (v- ip)(x) = 0. There- 
fore, with H3.31fl and (|3.32|) . we deduce that for m large enough : 

sup {v — ip m ) < = (v* - tp m )(T, x) < max {v* — <p m ). 

[s m ,T]xdB v (x) [s m .T]xdB n (x) 

In view of Lemma 13. II this proves that for m large enough, 

[s m ,T]xB,(i) is not a subset of M.(<p m ). (3.33) 
Step 3. From Ij3.28|l . notice that for all (t, x) £ [s m ,T] x B v (x), we have : 

for m large enough. In view of 13.33J1 and by definition of M(<p m ), we then may 
find some element (t,x) € [s m ,T] x B n (x) s.t. 

G(t,x,D x <p m (t,x),Dl<p m (t,x)) < 0. 

This is in contradiction with l|3.27|l . □ 
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3. 3. Some applications in finance 

3.3.1. Smooth-fit property of one- dimensional convex singular problem 

We consider the singular control problem introduced in Remark 12. 21 with linear 
dynamics : 

dX t = 8X t dt + "/X t dt + a t dt, (3.34) 

where 8 and 7 are constants with 7 > 0, and a is the control valued in R + . We 
consider the infinite horizon problem : 



v(x) — sup E 



(/(*?) 



Xctt ) dt 



where Xf is the solution to i|3.34|) starting from x at time 0. It is convenient to 
redefine control in terms L t = J a s ds, and in fact to enlarge the set of control 
processes to the set of nondecreasing, cadlag adapted processes L. This allows 
displacement in time of the control and also ensures the existence of an opti- 
mal control. This singular control problem models the irreversible investment 
problem for a firm. The process X represents the production capacity of a firm, 
which may increase it by transferring capital from another sector activity. 8 < 
is the depreciation rate of the capacity production, 7 its volatility, L t is the 
cumulated number of capital received by the firm until time t for a cost L XdL t , 
with A > 0, interpreted as a conversion factor from an activity sector to another 
one. / is the running profit function of the firm, assumed to be concave and 
finite, i.e. of linear growth condition, on (0, 00). 

The Hamiltonian of this one-dimensional convex singular control problem is : 



H(x,p,M) 



8xp + \ 1 2 x 2 M + f(x) if\-p>0 
00 if A — p < 0. 



Hence, the associated HJB variational inequality takes the form 
pv(x) — f3xv'(x) — -j 2 x 2 v / ' (x) — f(x) , A — v'(x) 



0. (3.35) 



One-dimensional convex singular control problems were largely sudied in the 
1980's and in particular their connection with optimal stopping problems. For- 
mally, the derivative of the value function for a singular control problem is 
the value function of an optimal stopping problem. In view of the smooth-fit 
principle for optimal stopping problems, which states that the associated value 
function should be differentiable, it is expected that the value function for a 
singular control problem should be twice differentiable. We show here how to 
derive simply the twice continuous differentiability for the value function of the 
above one-dimensional convex singular problem by using a viscosity solutions 
argument. 
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First, we easily check that by the linearity of the dynamics of the controlled 
process X and the concavity and linear growth condition of /, the value function 
v is concave with linear growth condition on (0, oo). Hence, v is continuous on 
(0,oo), and admits a left and right derivative v'_(x) and v' + (x) for each x > 0, 
with v' + {x) < v'_(x). We also see that 



v'_(x) < A, x > 0. 



(3.36) 



Indeed, fix a; > 0, and let us consider for any arbitrary control L and positive 
number < I < x, the control L consisting of an immediate jump of L of size 
/ at time 0. Then, the associated state process X satisfies : X x ~ l — X x and by 
definition of the value function, we have : 



v(x-l) > E 



-pi 



-pi 



(Kxr l )dt 

(f{X?)dt- 



XdL, 



XdL, 



XI. 



From the arbitrariness of L, this implies v(x — I) > v(x) — XI and so the required 



result 1|3.36[) by sending I to zero. By setting Xb = inf{x > 
using the concavity of v, we also easily see that 

AfT := {x > :v^{x)<X} = (x b ,oo), 

v is diffcrentiable on (0,Xb) with 

v'(x) — A, on (O^Xf,)- 



c) < A}, and 



(3.37) 



(3.38) 



We can also easily check that < xi, < 00. 

According to Theorem 13. II we know that v is viscosity solution to the HJB 
variational inequality H3.35JI . By exploiting furthermore the concavity of v in 
this one-dimensional problem, we show that v is actually a classical smooth 
solution. 

Theorem 3.3 The value function v is a classical solution in C 2 ((0, 00)) of 

Proof. Step 1. We first prove that v in C 1 on (0, 00). We argue by contradiction 
and suppose then on the contrary that v' + (x) < v'_(x) for some x > 0. Let us 
fix p e (v' + (x), v'_(x)) and consider the smooth function 



v(x) + p(x — x) (x — x) 2 , 

2e 



with e > 0. Observe that x is a local maximum of (v — ip £ ) with t^ e (x) = 
v (x). Since f' £ (x) = p < A by l|3.36|l and f' e '{x) = 1/e, the viscosity subsolution 
property implies : 



Pf{x) — [3xp H 7 2 x 2 

2e 



f{x) < 0. 



(3.39) 
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For e sufficiently small, we get the required contradiction, and so v',(x) = v'_(x). 

Step 2. By (|3.38|l . v is C 2 on (0,x&) and satisfies v'{x) — A, x £ (0,Xf,). From 
Step 1, we have NT = (x b , oo) — {x > : v'(x) < A}. Let us check that v is a 
viscosity solution of : 

pv(x) - (3xv'(x) - ]p 2 x 2 v"(x) - f(x) = 0, x<E(x b ,oo). (3.40) 

Let x G (xf,, oo) and ip a smooth test function C 2 on (x;,,oo) s.t. x is a local 
maximum of u — 99, with (v — <p)(x) — 0. Since f'(x) = v'(x) < A, the viscosity 
subsolution property of v to (|3.35(l implies : 

W (i)-^'(i)-i 7 2 iV(5)-/>) < 0. 

This shows that v is a viscosity subsolution of (|3.40() on (xi,,oo). The proof of 
supersolution viscosity of (|3.4()|l is similar. Consider now arbitrary x\ < X2 in 
(xb, 00), and the Dirichlet problem : 

pV{x) -PxV'(x) - ]p 2 x 2 V"{x) - f(x) = 0, xe( Xl ,x 2 ) (3.41) 
V(x 1 ) = v(x 1 ), V(x 2 )=v(x 2 ). (3.42) 

Classical results provide the existence and uniqueness of a smooth solution V, C 2 
on (xi, x 2 ), to (|3.41|l - (|3.42|l . In particular, this smooth solution V is a viscosity 
solution to l|3.4()(l on (xi,X2)- By standard comparison principle for viscosity 
solutions (here for linear PDE in bounded domain), we deduce that v = V on 
(x\, £2). From the arbitrariness of x\ and x 2 , this proves that v is C 2 on (x;,, 00) 
and satisfies (|3.40(l in the classical sense. 

Step 3. Recall that v is C 2 on (0, Xb) with v'(x) = A. It then remains to prove 
the C 2 condition of v in x\,. Let x 6 (0,Xb). The viscosity subsolution property 
H3.35JI of v in x applied with the smooth test function ip = v on (0, Xb), where 
if'(x) — A, (p"(x) = 0, implies that v satisfies 

pv(x) — p\x — f(x) > 0, < x < Xb, 

in the classical sense. By sending x to Xb, we get : 

pv(xb) - pXx b - f(x b ) > 0. (3.43) 

On the other hand, by the C 1 condition of v in Xb, we have by sending x to Xb 
in (|3~4Tty : 

pv(x b ) - (3\x b - f(x b ) = \l 2 x 2 b v"(x+). (3.44) 

From the concavity off, the r.h.s. of (|3.44|l is nonpositive, which combined with 
(|3.43|l . shows that v"{x^) = 0. This proves that v is C 2 in Xb with 

v"(x b ) = 0. 
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□ 

Once we have the C 2 regularity of the value function, we complete its explicit 
characterization as follows. We solve the ode ()3.38|1 on (0, Xb) and the ode H3.40J) 
in (xb, oo). We have then four unknown parameters : one coming from the first 
order ode (|3.38|> . two from the second order ode l|3.40fl . and the parameter Xb- 
These four parameters are explicitly determined by the linear growth condition 
of v and the C 2 condition of v in Xb- Moreover, the optimal state process and 
control are given by the reflected process X at the lower level Xb and by the 
local time L &t Xb- We refer to @HI for the details in a more general framework. 

The concavity property of v, resulting from the linear dynamics of X, is 
crucial for the C regularity of the value function : there are examples where 
the value function is not convex and not twice continuously differentiable. Some 
extensions to the two-dimensional case are studied in [HI] for a singular control 
problem arising from a transaction costs model, and the value function is shown 
to be twice differentiable under some assumptions. Singular control problems 
in multidimensional transaction costs model were recently studied in |50j but 
no regularity results are stated : only a viscosity characterization is given. It is 
an open question to know whether regularity results can be extended in such a 
multidimensional context. 

3.3.2. Superreplication in uncertain volatility model 

We consider the controlled diffusion : 

dX s = a s X s dW s , 

valued in (0, oo) (for an initial condition x > 0), and where the control a is 
valued in A = [a, a], < a < a < oo. To avoid trivial cases, we assume a > 
and a ^ oo. In finance, a represents the uncertain volatility process of the asset 
price X. Given a continuous function g with linear growth condition on (0, oo), 
representing the payoff function of an European option, we are interested in the 
calculation of its superreplication price given by, see e.g. |34j : 

v(t,x) = sup E [g(X%?)] , (t, x) e [0, T] x (0, oo) (3.45) 

Since the process X is a nonnegative supermartingalc, it is easy to see that v 
inherits from g the linear growth condition. In particular, v is locally bounded. 
The Hamiltonian of this stochastic control problem is 

H(x, M) = sup i ]-a 2 x 2 M \ , (x, M) e (0, oo) x R. 

ae[a,a] I 2 J 

We shall then distinguish two cases according to the finiteness of the volatility 
upper bound a. 

Bounded volatility 
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We suppose 

a < oo. 

In this regular case, the Hamiltonian H is finite on the whole domain (0, oo) x M, 
and is explicitly given by : 

H(x,M) = ^a 2 (M)x 2 M, 

with 



a(M) = 




if M > 
if M < 0. 



According to the results in paragraph s . 21 ( Thcorcms l3 .1113.21 and their follow- 
ing remarks), we have then the characterization result on the superreplication 
price v : 

Theorem 3.4 Suppose a < oo. Then v is continuous on [0,T] x (0, oo) and is 
the unique viscosity solution with linear growth condition of the so-called Black- 
Scholes-Barenblatt equation : 

- % - 1 &2 = °' M G [°' r ) x ( o ' oo )> ( 3 - 46 ) 

satisfying the terminal condition : 

v{T,x) = g(x), xe(0,oo). (3.47) 

Remark 3.6 When a > 0, there is existence and uniqueness of a smooth so- 
lution to the Black-Scholes-Barenblatt PDE (|3.46(l together with the Cauchy 
condition Ij3.47|l . and so v is this smooth solution. 



Unbounded volatility 

We suppose 

a = oo. 

In this singular case, the Hamiltonian is given by : 



H(x,M) = 



\a 2 x 2 M if - M > 



oo if - M < 0. 
According to Theorem l3.ll v is then a viscosity solution of 



dv 1 2 2 d 2v d 2y \ 

~di ~ 2- x 'dx 2 ' 



0, (t,x) G [0,T) x (0,oo).(3. 



Moreover, from Theorem 13.21 and its following remark, the terminal condition 
associated to the variational inequality H3.48J) is given by : 

v{T~,x) = g(x), xe(0,oo), 
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where g is the smallest function above g and satisfying in the viscosity sense 

d 2 q 

-q^( x ) ^ 011 (°>°°)- ( 3 - 49 ) 

If g were smooth, the previous inequality (|3.49|) characterizes the concavity of 
g. Actually, this is still true when l|3.49fl holds only in the viscosity sense, see 
|24| . Therefore, g is the upper concave envelope of g. We can then explicitly 
characterize the superreplication price v. 

Theorem 3.5 Suppose a — oo. Then v — w on [0, T) x (0, oo) where w is the 
Black-Scholes price for the payoff function g(x) : 

w(t,x) = E g(x^ , V(t,x) G [0,T] x (0,co), (3.50) 

in a Black-Scholes model with lower volatility a, i.e. {X l s ' x ,t < s < T} is the 
solution to 

dX s = g,X s dW s , t<s<T, X t = x. 
Proof. First, observe that w is the solution to the Black-Scholes PDE : 

~ '^~2- x dx I = °' on [0 ' T)x (0 ' w) ' (3 ' 51) 

with the terminal condition : 

w(T~,x) = w(T,x) = g(x). 

Since g is concave and it is well-known that Black-Scholes price inherits concav- 
ity from its payoff, we deduce that 

<9 2 u> 

-(t,x) > on [0,T) x (0,oo). 



dx 2 

Together with the previous equality (|3.51(l . this proves that w is solution to 

, -— \ = 0, (t,x) G 0,T) x (0,oo). 



dt 2- dx 2 ' dx 2 , 

Finally, since w and w satisfy the same terminal condition, we conclude with a 
comparison principle on the PDE (|3.48ll . □ 

Other variations and extensions to financial models with portfolio constraints 
or transaction costs are studied in |23] and f° r the calculation of superrepli- 
cation price by a viscosity solutions approach. 
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4. From stochastic control to backward stochastic differential 
equations 

As already mentioned in paragraph 12.31 there is a strong relation between the 
stochastic maximum principle and backward stochastic differential equations, 
BSDE in short. A BSDE is usually written in the form 

dY t = -f{t,Y t ,Z t )dt + Z t dW t , Y T = £. (4.1) 

Here F = (J-t)o<t<T is the filtration generated by W, the terminal condition 
£ is an ^V-measurable random variable, and a solution to (|4.1|l is a pair of 
F-adapted processes (Y, Z) s.t. 

Y t = C + J t f(s,Y s ,Z s )ds- J Z s dW s , 0<t<T. 

Existence and uniqueness of a solution (Y, Z) are usually stated under a Lip- 
schitz condition on /, see Pardoux and Peng [73] • Some extensions to the 
quadratic growth case are studied in and more recently in 02]- We re- 
fer to [S3 for an account on the theory of BSDE and to jSH] for its applications 
in finance. 



4-1- Stochastic maximum principle and BSDE 

We state a sufficient maximum principle for a stochastic control problem in the 
framework described in Section 0: a controlled state process 

dX s = b(X s ,a s )ds + a(X s ,a s )dW s , X a = x, 

and a gain functional to maximize : 



J (a) = E 



f(t,X t ,a t )dt + g(X T ) 



Moreover, the coefficients b(t,x,u>), a(t,x, u>) of the state process as well as the 
gain functions f(t,x,u>) and g(x,u>) may depend on u> : b(t,x,cj), a(t,x,uj) 
and f(t,x,uj) are ^-measurable for all (t, x) and g(x, u) is JF T -measurable. For 
simplicity of notation, we omit the dependence in w. We consider the generalized 



pnxd 



Hamiltonian H : [0, T] x M" x A x W 

H(t,x,a,y, z) — b(x, a).y + tr(a'(x, a)z) + /(i, x, a). 



(4.2) 



We also omit the dependence of H in u. Notice the difference with the Hamilto- 
nian (|2.7J) introduced in the dynamic programming approach. We suppose that 
H is differentiable in x with derivative denoted D X TL, and for each a e A, we 
consider the BSDE, also called the adjoint equation : 



dY t = -D x H(t,Xt,a t ,Y t ,Z t )dt + Z t dW t , Y T - D x g{X T ). (4.3) 



The following sufficiency theorem is proved e.g. in |89| . 
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Theorem 4.1 Assume that g is concave in x. Let a € A and X the associated 
controlled state process. Suppose there exists a solution (Y, Z) to the correspond- 
ing BSDE <|4.3[) such that almost surely : 

H(t,X t ,a t ,Yt,Zt) = m&xH(t,X t ,a,Yt,Zt), < t < T, (4.4) 

and 

(x,a) — > TC(t, x, a, Y t , Z t ) is a concave function , (4-5) 
for all t £ [0, T]. Then a is an optimal control, i.e. 

J(a) — sup J{a). 

4-2. Application : linear- quadratic control problem and 
mean-variance hedging 

Mean-variance hedging is a popular criterion in finance : this is the problem 
of finding for a given option payoff the best approximation by means of a self- 
financed wealth processes, and the optimality criterion is the expected square 
error. This problem has been solved in high generality for continuous semi- 
martingale price process and general filtration by martingales and projection 
techniques, see (ZH) aud [EH! for recent overviews. On the other hand, in a series 
of recent papers, this problem has been reformulated and treated as a linear- 
quadratic (LQ) control problem at increasing levels of generality, mostly for Ito 
processes and Brownian filtration, see e.g. jSU, [531, EH- 52- When the coef- 
ficients are random, which generalize the classical case of deterministic linear- 
quadratic control problem, the adjoint equations lead to a system of BSDEs. 
This gives some new insight and results on general LQ control problem, and 
provides also fairly explicit results for the mean-variance hedging. 

The problem is formulated as follows. For simplicity of notation, we consider 
the one-dimensional case. The dynamics of the linear controlled state process 
is : 

dX t = (r t X t + a t [it)dt + crtCttdWt, X = x, 

and the objective is to minimize over controls a valued in R, for a given square 
integrable ^r-measurable random variable £, the quadratic functional : 

J (a) = E[X T -tf. 

The coefficients r t , fi t and at are assumed to be bounded and adapted with 
respect to F, the natural filtration of W. We also assume a uniform ellipticity 
condition on a : at > £ a.s. for some e > 0. In finance, X is the wealth process 
of a self-financed portfolio a, representing the amount invested in a stock asset 
of excess return fi t with respect to the interest rate , and volatility er t . £ is the 
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payoff at time T of an option that we are trying to approximate by X according 
to the quadratic error criterion. 

We show how to solve this optimization problem by applying the stochastic 
minimum principle in Theorem 14. II Observe that the generalized Hamiltonian 
for this stochastic LQ problem takes the form 

H(x,a,y,z) = {r t x + afj. t )y + cr t az, 

and the adjoint BSDE 14.3J1 is written, given a G A : 

dY t = -r t Y t dt + Z t dW t , Y T = 2(X T -0- ( 4 -6) 

Suppose a G A is a candidate for an optimal control, and denote X, (Y, Z) 
the corresponding processes. Since the Hamiltonian is linear in a, we see that 
conditions (|4.4() and (|4.5() will be satisfied iff : 



H t Y t + (j t Z t 



0, < t < T, a.s. 



(4.7) 



Due to the linear-quadratic formulation of the problem, we are looking for a 
solution (Y, Z) to Ij4.6jl satisfying Ij4.7jl in the form 



Y 



P t (X t -Q t ), 



(4.8) 



where P and Q are adapted processes satisfying the terminal condition : Pp — 
2 and Qt = Searching for Ito processes and proceeding by identification, we 
see after straightforward calculation that if there exist pairs (P, A) and (Q, T) 
solutions to the BSDEs : 



dP t = -Pt. 



at 



Pt 



dt + K t dWt 



dQ t 



r t Qt + -r t 



T t dW t , Q T = t, 



(4.9) 



(4.10) 



then (Y, Z) is solution to the BSDE (|4.6(l with the mimimum condition l|4.7|l . 
Existence and uniqueness of a solution to H4.9|I - (|4.1U|I are proved in More- 
over, the optimal control is given by : 



a i 



1 

0t 



at 



A* 
Pt 



(Xt 



r* 



(4.11) 



Finally, the value function is simply derived by applying Ito's formula to P t (X t - 
Qt) 2 /2 for a G A : 



J (a) = E[X T -t] 2 = E 



P T (X T - Q T f 



-Po(x-Qo) 



2 



/ P t (J 2 {at - atfdt 
Jo 



This proves again that the optimal control is indeed given by (|4.11|) and the 
value function is equal to : 

1 



J(a) 



-Pn(x 



1 
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4-3. Forward /backward stochastic differential equations and control 

Suppose that the conditions of the sufficiency theorem 14. f I are satisfied and 
denote by X the corresponding optimal controlled state process, and (Y, Z) the 
associated adjoint dual variables. By considering the function 

H(t,x,y,z) = sup TC(t, x, a, y, z) — sup [b(x, a).y + tr(<j (x, a)z) + f(t, x, a)) , 

we then see, by the envelope theorem, that (X, Y, Z) satisfy : 

dX t = H y {t,X u Y t ,Z t )dt + H z (t,X u Y u Z t )dW u 
dY t = -H x (t,X u Y u Z t )dt + Z t dW u Y T = D x g{X T ). 

where (H x H y H z ) stands for the gradient of H with respect to the (x y z) vari- 
ables. This is a special case of coupled forward-backward stochastic differential 
equation (FBSDE) written in the form : 

dX t = b(t,X u Y u Z t )dt + <j(t,X t ,Y t ,Z t )dW t , (4.12) 
dY t = -f(t,X t ,Y u Z t )dt + Z t dW t , Y T = h{X T ). (4.13) 

A solution is a triple (X, Y, Z) of adapted processes with X a solution to the 
forward s.d.e. I|4.12|l and (Y, Z) a solution to the bsde (|4.13|) . Coupled FBSDEs 
extend BSDEs since the backward component (Y, Z) appear in the coefficients 
of the forward component X. The study of FBSDEs generates an important 
stream of research. There are essentially two methods for proving the solvability 
of BSDE : pure probabilistic methods like the method of continuation, see [2], 
|4"5] , and combined methods of PDE and probability like the four step scheme 
developed in |65| . or more recently extensions in |74| . |27j. The connection with 
PDE is an extension of the Feynman-Kac formula : it states that if there exists a 
smooth solution with suitable Lipschitz and growth conditions to the quasilinear 
PDE 

~~&t ~ Xl Vl X ' U ' DxV ^' DxV ~ 2 tr ( crc7/ (*' Xl Vl Xl u ' DxV )) D x v ) 
-f(t,x,v,z(t,x,v,D x v)) = 0, (t, x) e [0, T) x R n (4.14) 
v(T,x) = h(x), x e M™, (4.15) 

where z(t,x,y,p) satisfies z(t,x,y,p) = p' o(t,x,y, z(t,x,y,p)), then the triple 
(X, Y, Z) determined by 

dX t = b(t,X t )dt + a(t,X t )dWt, 

where 

b(t,x) — b(t,x,v(t,x),z(t,x,v(t,x),D x v(t,x)) 

5{t,x) = a(t,x,v(t,x),z(t,x,v(t,x),D x v(t,x)), 
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and 

Y t = v(t,X t ), Z t = z(t,X t ,v(t,X t ),D x v(t,x t )) 

is a solution to the FBSDE H4.12fl ~ H4.13fl . Conversely, if there exists a solution 
(X t < x ,Y t ' x ,Z t < x ) solution to the FBSDE ^T^ - ftT^ with X*' x = x, then the 
function v(t,x) = Y^' x is a viscosity solution to If4.14fl - lf4.15fl . We refer to the 
lectures notes j^El for an account on the theory of FBSDE and their applications. 

5. Numerical issues 

Explicit analytic solutions to stochastic control problems are rare and one has 
to resort to numerical approximations. This has led to an important area of 
research on numerical methods for stochastic optimization. Motivated in partic- 
ular by problems in quantitative finance, one faces challenging numerical prob- 
lems arising in portfolio optimization in high dimension, under partial informa- 
tion, under transaction costs ... We distinguish two types of numerical methods 
for stochastic control problems : purely deterministic and probabilistic methods. 
We briefly survey some of the developments and advances in these fields. 

5.1. Deterministic methods 

Purely deterministic methods from numerical analysis consist in discretizing 
the nonlinear HJB PDE satisfied by the value function of the stochastic con- 
trol problem. The discretization by finite difference methods or finite elements 
methods leads to an approximation of the value function at the points of the 
space-time grid. Convergence of the numerical scheme is proved rigorously by 
stability results for viscosity solutions, see |7]. There are also recent results on 
convergence rate for the finite difference approximation of the HJB equation in 
|K7] and [Hj- Computational methods are studied e.g. in p. For some illustrations 
in financial problems, see e.g. or [HHj- From a computational viewpoint, the 
limitation of purely deterministic methods is the dimension of the state variable 
and in practice, calculations are done for low dimensions, say 1 or 2. 

5.2. Probabilistic methods 

Probabilistic numerical methods for optimization problems considerably devel- 
oped over these years. A classical method, based on the dynamic programming 
principle, is the Markov chain approximation introduced by Kushner |58| . see 
also the book for an account of this method. This consists in approxi- 
mating the original continuous-time controlled process by a suitable controlled 
Markov chain on a lattice satisfying a local consistency condition. The stochas- 
tic control problem is then numerically solved from the dynamic programming 
principle applied to the Markov chain approximation, which leads to a back- 
ward recursive algorihm. The finite difference scheme is a typical example of 
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a numerical scheme for an approximating Markov chain with nearest neighbor 
transitions. In this method, the required stability condition may be very restric- 
tive in the case of control-dependent diffusion. On the other hand, the lattice 
is regardless of the structure of the Markov chain. Moreover, its size is grow- 
ing exponentially with the dimension. So, although the Markov chain is easily 
implemented, calculations can be done in practice only for low dimensions as 
in the purely deterministic approach. To overcome the dimensionality problem, 
Pages, Pham and Printems [7J propose a probabilistic method based on opti- 
mal quantization methods for numerically solving stochastic control problems 
in dimensions larger than 3. Like in usual probabilistic methods, we start from a 
time discretization of the controlled problem : the process (X t ) is approximated 
by its Euler scheme denoted Xk at time tk = kT/n. Then, and in the spirit of 
the Markov chain approximation method, the Euler scheme is approximated at 
every date A: by a process Xk, taking finitely many states. The optimal quantiza- 
tion approach consists in finding the best approximation according to L p -norm 
of Xk by Xk ■ This is achieved by stochastic gradient descent method based on 
Monte-Carlo simulations of Xk, which also allows to estimate the probability 
transitions of the grid points. The control problem is then numerically imple- 
mented from the backward dynamic programming formula. The main interest 
in optimal quantization is that given a total number of points to be dispatched 
among all grids, one gets optimal grids with respect to the distribution structure 
of the original process. This method is numerically tested in dimension 3 for the 
mean-variance hedging problem with stochastic volatility. We refer to [72] for 
an account of other applications of quantization methods to numerical problems 
in finance. 

Another approach is to solve numerically the FBSDE associated via the max- 
imum principle to the stochastic control problem. Several approximation meth- 
ods to FBSDEs were proposed in the literature. We mention for BSDE the 
Markov chain approximation in [53] . regression methods in 03], quantization 
methods in |3J, the Monte-Carlo Malliavin method in [TH] , and for FBSDE the 
four-step scheme in |2H], and quantization method in 

6. Conclusion 

Many theoretical and numerical advances have been recently realized in the 
field of stochastic control. They also contribute to the fields of nonlinear partial 
differential equations and backward stochastic differential equations. Present 
developments include modelling with jump-diffusion processes, which lead to 
integrodifferential equations, and FBSDEs with reflecting barriers arising typ- 
ically in optimal stopping problems. From a numerical viewpoint, challenging 
problems are the search for fast and efficient schemes for stochastic control 
problem in high dimension arising in quantitative finance. In this direction, par- 
tial observation problems, which lead to infinite-dimensional problem, are very 
few studied numerically and represent an important stream of research. Clas- 
sical control problems in finance are optimal portfolio allocation but there is a 
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large potential applicability for other economics and finance problems such as 
in contract or game theory, credit risk or liquidity risk models. In return, these 
questions should raise new developments in terms of mathematical theories, still 
in accordance with meaningful applied problems. 
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